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We have studied terahertz (THz) emission from arsenic-ion implanted GaAs both experimentally 
and using a three-dimensional carrier dynamics simulation. A uniform density of vacancies was 
formed over the optical absorption depth of bulk GaAs samples by performing multi-energy im- 
plantations of arsenic ions (1 and 2.4 MeV) and subsequent thermal annealing. In a series of THz 
emission experiments the frequency of peak THz power was found to increase significantly from 
1.4 to 2.2 THz when the ion implantation dose was increased from 10^"^ to 10^^ cm~^. We used a 
semi-classical Monte-Carlo simulation of ultra-fast carrier dynamics to reproduce and explain these 
results. The effect of the ion-induced damage was included in the simulation by considering carrier 
scattering at neutral and charged impurities, as well as carrier trapping at defect sites. Higher va- 
cancy concentrations and shorter carrier trapping times both contributed to shorter simulated THz 
pulses, the latter being more important over experimentally realistic parameter ranges. 

PACS numbers: 78.47.+p, 78.20.Bh, 61.72.Vv, 68.55.Ln DOI: 10.1103/PhysRevB. 70.235330 



I. INTRODUCTION 

The technique of terahertz time-domain spectroscopy 
(THz-TDS) has seen wide application in recent years. By 
resolving carrier dynamics on ultra-short time scales THz- 
TDS has proven invaluable in areas of physics as disparate 
as investigating the energy gap in high Tc superconductors-, 
determining low-energy vibrational modes in oligomers-, and 
observing the onset of quasi-particle screening in GaAs? . 

Such studies have become possible through the develop- 
ment of large bandwidth, high power THz emission sources 
based on the ultra-fast separation of photoexcited carriers 
in semiconductors. The resulting change in current pro- 
duces an electromagnetic transient, typically a single or half- 
cycle pulse, with a broad spectrum covering THz frequencies. 
Charge separation can be the result of differing electron and 
hole mobilities (the photo-Dember effect), or alternatively 
may result from an electric field accelerating electrons and 
holes in opposite directions. In the later case, the electric 
field can be internal, such as the surface depletion field near a 
semiconductor surface, or externally applied as in a photocon- 
ductive switch (PCS). The study of surface emitters can help 
to elucidate aspects of ultra-fast carrier dynamics, such as 
the enhancement of THz emission under a magnetic field^*^, 
and the relative importance of the photo-Dember and sur- 
face field charge separation mechanisms in GaAs and InAa-i-. 
Additionally, surface emitters can produce the well-coUimated 
THz beams^ desirable for many THz imaging applications. 

In the effort to improve THz emitters, one method for 
achieving a greater bandwidth is to shorten the electric field 
rise time, normally limited by the duration of the exciting 
infrared pulse. Another approach is to reduce the duration 
of the electric field transient after excitation by choosing a 
semiconductor with a short carrier lifetime, and/or a high 
momentum scattering rate. One suitable material is low- 
temperature-grown (or LT-) GaAs, in which As precipitates 



and trapping centres result in a short (< 1 ps) carrier life- 
time, enabling THz emission extending to over 30THz^. Al- 
ternatively the technique of ion implantation can be used, in 
which incident energetic ions damage the crystal structure of a 
semiconductor. The principal defect in arsenic ion-implanted 
GaAs (GaAsiAs""") is thought to be the antisite donor de- 
fect AsGa, where As replaces Ga on some Ga lattice sites. 
Deep-level transient spectroscopy has found that the defect 
energy levels in GaAs:As^ lay ~0.3eV below the conduction 
bandpi-. Such deep-level defects can decrease THz pulse dura- 
tion (increase THz bandwidth) by acting as carrier trapping 
and scattering centres, and result in carrier lifetimes of the 
order of 0.1 psiiii^. 

The controllable ion dosage permits ion-implanted GaAs 
samples with reproducible properties, in contrast to LT-GaAs 
where the difficulty of temperature control results in samples 
with nominally the same growth parameters having varying 
properties. Ion-implantation has the further advantage of 
generating depth-dependent damage profiles (by choosing the 
ion type, dose, and energy) in specific sample areas. By per- 
forming a post-implant anneal the resistivity of GaAsiAs""" 
can be increased to levels comparable to LT-GaAs—*—, as 
required for photoconductive switches. 

Both surface field— and photoconductive antennai^ ion- 
implanted THz emitters have been previously studied, and 
were found to produce THz radiation at slightly higher fre- 
quencies than semi-insulating GaAs^^. Additionally, ion- 
implanted photoconductive antennae have been used as 
detectors^'^, but no bandwidth improvement could be seen. 
These studies utilised low energy ions, typically 200 keV, re- 
sulting in implantation depths of only ~Q.l/im. However, as 
the infrared absorption depth is ^^0.8 pm in GaAs only ^^10% 
of photocarriers are generated within this distance, and the 
non-implanted layer dominates THz emission. 

In this paper we investigate GaAs implanted with dual, 
high energy (IMeV and 2.4 MeV) As ion doses, creating an 
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FIG. 1: (Colour online) Depth dependence of vacancies (in- 
cluding recoil vacancies) calculated using SRIM software" for 
the sample with highest ion dose. The 1 MeV, 2.5 x 10^^ cm~^ 
dose (dash-dotted line) and 2.4 MeV, 1.0 x lO" cm"^ dose 
(dashed line) create an approximately uniform damage profile 
(solid line) extending over 1 /xm from the surface. By scaling 
the concentration according to the dose, the distributions for 
lower dose samples can be obtained. 



approximately uniform damage profile extending ~l/im into 
the semiconductor. In Section [H] we describe the terahertz 
emission from samples irradiated with various ion doses, and 
discuss the effect of a post-implant thermal annealing step. 
Simulation results of THz emission from ion-implanted GaAs 
obtained from a three-dimensional carrier dynamics model'' 
are presented in Section IITTI 



II. EXPERIMENT 

The Australian National University 1.7 MV tandem accel- 
erator was used to irradiate GaAs samples from the same 
ingot with As""" ions. By employing a two-stage implant 
of 1 MeV and 2.4 MeV (implanted at room temperature) an 
approximately uniform damage profile >lfim in depth was 
created (Fig. in 7 samples at incident ion doses ranging 
from 2.5 x 10^^ cm'^ to 2.5 x 10^^ cm"^ (1 MeV implant) and 
1 x 10^^ cm"^ to 1 X lO^^cm-2 (2.4MeV implant). A du- 
plicate set of samples was produced that were annealed for 
30min. at 500° C under AsHs. Calculations using the Stop- 
ping Range of Ions in Matter (SRIM) softwarei^i predict that 
the implanted ion concentration in the lowest dose sample 
is ~ 1 X lO^'^cm"^. Assuming that 99% of Asca antisites 
are instantly annealed^^ due to the implantation being per- 
formed at room temperature, we estimate an n-type doping 
concentration of ~ 1 x 10^^ cm~'^ for the lowest dose sam- 
ple. Damage distributions were also calculated using SRIM, 
and verified using X-ray diffraction. Approximately 10,000 
vacancies were found to be produced per As^ ion. 

The surface field THz emission from each GaAs:As+ sample 
was studied using THz-TDS. Infrared pulses from a mode- 
locked 10 fs Ti-sapphire laser (A — 800 nm, AA = 135 nm, 
75MHz repetition rate, ~ 400 mW typical beam power) were 
incident on the sample at 45°. A lens focused the infrared 
pulses to a Gaussian beam- waist of cr = 90 fim. The sam- 
ple was positioned ~lcm from the pump beam focus, with 
a Gaussian beam-waist of cr = 160 ^im chosen to provide a 



coUimated THz beam^. The THz radiation was collected by 
an off-axis parabolic mirror placed at its focal length from the 
pump beam focus, and refocused by an identical mirror onto 
a 0.2 mm thick (110) ZnTe crystal. A delayed probe beam 
was split from the pump using a 10% beam-splitter, and also 
illuminated the ZnTe. The change in the probe's polarization 
was measured as a function of delay using an electro-optic 
sampling (EOS) system consisting of a WoUaston prism, A/4 
waveplate and balanced photodiodes. The pump beam was 
chopped at ~ 1 kHz, enabling the voltage across the detec- 
tor to be measured by a lock-in amplifier (200ms time con- 
stant), from which the THz electric field was calculated—. 
Following the procedure of Zhao et. al^ we determine that 
our EOS system operates in the shot-noise limit, with a 
noise-equivalent THz electric field -Enef — 0.01 Vm^^ consis- 
tent with the experimentally observed root-mean-square noise 
EuMS ^ 0.01 0.02 Vm~\ 

The THz electric field emitted by each of the 7 annealed 
samples was thus measured, with signal-to-noise ratios of 
above 150:1. The maximum THz electric field -Ethz was 
of the order of 3Vm~^, less than we observe for our unim- 
planted GaAs reference sample (4 Vm~^), and an InAs sample 
(216 Vm~^) with the same experimental setup. For compar- 
ison, the maximum THz electric field strength reported in 
the literature from a photoconductive antenna is 9500 Vm~^ 
(Ref.E). 

^-THz was not found to vary systematically with 
ion dose, perhaps because of variations in the surface poten- 
tial of the samples, which can greatly alter the surface field 
strength. 

Typical THz electric fields are shown in Fig. |5| for low 
(2.5 x 10^^ cm"^ at IMeV, 1 x lO" cm~^ at 2.4 MeV) and 
high (1.25 X 10^^ cm"^ at IMeV, 5 x 10^^ cm'^ at 2.4 MeV) 
implant doses. At higher implant doses we observe that the 
THz pulse duration is shorter, and that the negative peak af- 
ter the principal peak has both larger magnitude and smaller 
period (Fig. |5Ja)). In Fig. |5Jb) the Fourier transforms of the 
time-domain electric field show that at greater ion doses THz 
emission shifts to higher frequencies. 

The post-implant thermal annealing step (500°C, 30min.) 
was observed to result in larger peak THz fields, and greater 
power at low frequencies, as shown in Fig.j^for the lowest dose 
sample (rid = 2.5 x 10^^ cm"^ at 1 MeV, rid = 2.5 x 10^^ cm"^ 
at 2.4 MeV). Annealing removes defects, at least partly re- 
pairing the damaged crystal structure. Because defects scat- 
ter carriers, and can also trap conduction band electrons, we 
can qualitatively predict that annealing will increase the pulse 
width, decreasing the relative THz emission at high frequen- 
cies. 



III. SIMULATION 

In order to model the ultra-fast carrier dynamics leading 
to terahertz emission from semiconductor surfaces in such a 
way that reproduces the experimentally observed features it 
is appropriate to employ a three-dimensional Monte Carlo 
simulation, such as presented in Ref. 0. This model included 
the interactions between extrinsic and photogenerated carri- 
ers in r and L-valleys, plasmon and magnetoplasmon effects, 
and the dielectric-air interface, and has been used to distin- 
guish between the diffusion (photo-Dember) and surface field 
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FIG. 2: (Colour online) Measured THz electric field -Eth: 
emitted via the surface field mechanism from GaAsiAs"'" 
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pies implanted with low doses (solid lines, 1 x 10 cm at 
2.4 MeV and 2.5 x lO" cm"^ at 1 MeV) and high doses (dot- 
ted lines, 5 x 10^^ cm'^ at 2.4 MeV and 1.25 x 10^'^ cm'^ at 
IMeV). The time and frequency domain THz electric fields 
are shown in (a) and (b). The peak electric fields and trailing 
negative peaks are -1-3.0, — 1.5 Vm~^ and +3.7, — 2.0Vm~^ for 
the low and high dose samples respectively. 
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FIG. 3: (Colour online) THz emission from unannealed (dot- 
ted lines) and annealed GaAs:As+ (solid lines) for the lowest 
implant dose (1 x 10^^ cm"^ at 2.4 MeV and 2.5 x 10^^ cm"^ at 
1 MeV). Annealing causes an increase in maximum THz field 
from 1.6 Vm~^ (unannealed) to 3.0 Vm~^ (annealed), but a 
decrease in the frequency of peak power from 2.1 THz (unan- 
nealed) to 1.6 THz (annealed). 



mechanisms of charge separation in InAs and GaAs^*-, and 
to provide a quantitative explanation of the enhancement of 
THz emission under a magnetic field&. 

In this section we present simulations of the THz emission 
from GaAs:As^ using a version of the model of Ref. 'Qi, ex- 
tended to include the effects of the significant concentrations 
of vacancies and trapping defects. All simulations used the 
physical properties of GaAs as given in Table 1 of Ref.|^ and 
the laser parameters of Sec. ^ Defect concentrations were 
assumed to be uniform over the range of infrared absorption, 
as is approximately the case in Fig. Q Altering the charged 
impurity (donor) concentration was not found to significantly 
alter the spectral shape of THz emission. 



A. Neutral impurity scattering 

When As^ ions collide with atoms in GaAs both the in- 
cident ion and the recoiling target atoms introduce vacan- 
cies, which act as carrier scattering centres. Calculations us- 
ing SRIMi- suggest that approximately 10,000 vacancies are 
produced per incident ion, with the damage profile of FiglTl 
For our unannealed samples this corresponds to vacancy con- 
centrations ranging from A'vac ~ 10^^ cm~^ to ~10'^^cm~^. 
Electron paramagnetic resonance experiments have shown— 
that the concentration of arsenic antisite defects is about 
10^* cm~^ after a 500°C, 30min. anneal, for ion doses of 
lO^^cm"^. Since A'^vac for our annealed samples therefore 
varies from 10^^ cm~'^ to 10^* cm"'', carrier- vacancy scatter- 
ing cannot be disregarded as it normally is at lower defect 
concentrations. 

Vacancies are included in the simulation by assuming that 
carriers scatter elastically from an electrically neutral spher- 
ical square well of width R, as defined by the potential 

V — —Vo for carrier- vacancy distances r < R and V = for 
r > R. Following the partial- wave approach taken in Ref. I20I 
the momentum scattering rate Wvac can be calculated from 
the I = wave cross-section a using W^ac ~ a-N^i^cV, where 

V is the electron velocity before (and after) scattering. The 
resulting expression is: 



VKvac = 



47rfiiV„ 



(k cot kR — a cot aR)^ 



m*k (fc2 -I- a2 cot^ ai?)(l + cot^ fci?) 



where a is given by: 



2m'- {E + Vo) 



ft2 



(2) 



The final state carrier is scattered uniformly into 4n stera- 
dians, because of the spherical symmetry of the well. Taking 



the values of 7? = 3 nm and Vo 
tering rate ranging from ~ 10^' 



300 meV™ produces a scat- 
at iVvac = 1 X 10^^ cm~^ 



to 



10^ 



at iVvac = 1 X 10^ 



To enable a comparison between these momentum scat- 
tering rates and those from mechanisms previously included 
in the simulation, we plot in Fig. |1| the average momentum 
scattering rates of electrons in the F- valley as a function of 
time t after the infrared pulse for NdyN^ac = 1 x 10^^ cm~''. 
The scattering rates after f = Ops are averaged over all pho- 
toinjected and extrinsic electrons, whereas before the infrared 
pulse arrives only extrinsic electrons contribute. The time de- 
pendence of the scattering rates is a consequence of the evo- 
lution of the carrier energy and concentration distributions. 
The greatest contributions to the total rate are attributed to 
electron-hole scattering and charged impurities, which both 
decrease eit t — since they are inversely proportional to the 
local carrier density^. LO-phonon emission and absorption 
and TO-phonon emission and absorption followed by a jump 
into the L- valley initially contribute ^-^10% of the total rate, 
but become less significant at later times when the average 
carrier energy has decreased. Acoustic phonon and neutral 
vacancy scattering provide only ~1% of the total rate. 

While Fig. |1] provides some insight into the significance of 
the various momentum scattering mechanisms, the angular 
distribution of each mechanism must also be considered. At 
A'^vac = 1 X 10^^ cm~^ the vacancy scattering rate is ~ 6 x 
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FIG. 4: (Colour online) Momentum scattering rates for elec- 
trons in the F-valley, with donor and vacancy concentrations 
of 1 X lO^^cm""^ each. The rate for each scattering mech- 
anism is averaged over all particles in the simulation, and 
is plotted as a function of time after laser excitation. The 
largest contributions to the total rate (solid line) come from 
electron-hole scattering (circles) and charged impurity scat- 
tering (dots). LO-phonon emission (triangles) and absorption 
(dashed line) are also significant. Acoustic phonon scattering 
(squares) and vacancy scattering (diagonal crosses) produce 
lower rates. Least important are TO-phonon absorption (di- 
amonds) and emission (vertical crosses), where in both cases 
the electronic final state is in the L- valley. The infrared pulse 
arrives ai t = Ops, and produces a step in most scattering 
rates. 



FIG. 6: (Colour online) Simulated THz electric field from bulk 
GaAs at carrier trapping times of Tc = 10 ps and Tc = 0.1 ps. 
In (a) the unfiltered simulated Tc = 0.1 ps pulse (dashed line) 
has shorter duration than the 10 ps (solid line). Their Fourier 
transforms are shown in (b). 



B. Carrier-trapping defects 

The effect of trapping defects on THz emission was included 
in our model by introducing an exponential decay of the num- 
ber of photoexcited carriers n as a function of time t after the 
infrared pulse according to the equation; 

n{t) = n(0)e"'/"'^ (3) 
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FIG. 5: (Colour online) Simulated THz emission in (a) the 
time domain and (b) the frequency domain with no vacancies, 
and for vacancy concentrations of 1 x lO^^cm"^, 5 x lO^'^cm"^. 
The FWHM of the spectra increase with vacancy concentra- 
tion. Coherent phonon oscillations^''^ are not observed because 
the response of the crystal lattice is not included in the sim- 
ulation. 



10^^ s~^, i.e. only 10% of the total rate, yet as can be seen 
in Fig. EI a) the THz pulse duration is reduced. This is a 
consequence of the uniform angular distribution for vacancy 
scattering: the carrier direction is altered more significantly 
than in the other mechanisms. 



where Tc is the carrier trapping lifetime. Time-resolved pho- 
toluminescence experiments of GaAs:As^ implanted with a 
single 10^® cm~^ dose of 2 MeV ions have measured the car- 
rier trapping time to be as short as Tc = 0.1 ps— . It is im- 
portant to distinguish Tc from the carrier recombination time, 
which is ~4 pa—. We assume that once trapped, carriers can- 
not escape the defect via thermal excitation, and that trapped 
carriers do not alter free carrier states. When the number of 
carriers n is reduced on time scales of < 1 ps the electric field 
decay time after the pulse can be shortened, as shown in Fig. 
|n| resulting in a greater bandwidth. At Tc = 0.1 ps the FWHM 
from the simulation is 7.8 THz, nearly 80% larger than at a 
carrier trapping time of 100 ps (FWHM = 4.4 THz). 



IV. COMPARISON OF EXPERIMENT AND 
SIMULATION 

The spectrum of the simulated THz emission has a larger 
magnitude at both low and high frequencies than the exper- 
imental data (Fig. IT^b)). This can be accounted for by the 
effect of the EOS system and the parabolic mirrors used to 
collect the THz radiation. A frequency-domain picture of 
EOS with ZnTe has been presented by Gallot et. al. 
that models the effect of the detection apparatus using three 
complex filters, namely: the spectrum of the autocorrelation 
of infrared pulse (which has negligible effect), the frequency- 
dependent electro-optic susceptibility x'^' of ZnTe, and the 
mismatch between the infrared group velocity and THz phase 
velocity in ZnTe. The resulting low-pass filter causes the high 
frequency components to reduce in magnitude, as shown in 
Fig.lTJb). While the slow drop-off in simulated electric field at 
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FIG. 7; (Colour online) Comparison of simulated THz emis- 
sion from GaAs with n-type doping A'd = 1 x lO'^^ cm~"^ us- 
ing the low-dose sample data from Fig. |21 In (a) the simu- 
lated electric field (dots) is seen to have a sharp rise ~ 10 fs 
in duration, and to drop sharply after ~ 0.1 ps, producing 
the broad spectrum in (b). When low-pass and high-pass fil- 
ters are applied (taking account of the dispersion in ZnTe, 
and diffraction-limited collection respectively, as described in 
the text), the filtered spectrum (dashed line) can be seen 
to match the experimental spectrum (solid) extremely well. 
When transformed back into the time-domain, the filtered 
simulated electric field (dashed line of (a)) matches for times 
> 0.2 ps. 



increasing frequencies suggests that large bandwidths are pos- 
sible from surface field THz emitters, experimentally a thin 
(~ lO/xm) ZnTe crystal is necessary to observe such frequency 
components, which would result in a prohibitively small sig- 
nal. 

While the experiment was designed to collect THz radia- 
tion with high efficiency, low frequency THz radiation emit- 
ted from the sample surface will be diffracted beyond the 
collection capability of the parabolic mirrors. By treating 
the two parabolic mirrors as a single thin lens. Cote et. al. 
have calculated the high-pass filter corresponding to this 
effect^®. They found that the approximate filter function 
-fhigh('^) = erf^ (2a; i?o"/c/), where the THz radiation has a 
Gaussian beam-waist a, and is collected by two parabolic mir- 
rors each of focal length / and radius R. For our experimental 
geometry (/ = 15 cm, R = 2.5 cm, a — 0.16 mm) this filter 
has a steep rise from to 1 between and ~ 1.5 THz, which 
when applied to the simulated data causes a reduction in low 
frequency components of the electric field. 

Multiplying the simulated spectrum by these filters in the 
frequency domain produces a high accuracy match to the ex- 
perimental spectra (Fig. EJb)). Taking the inverse Fourier 
transform of the filtered simulated spectrum results in a time- 
domain trace that compares adequately with the experimen- 
tal time-domain data (Fig. [7[a)). The low-pass (ZnTe) filter 
causes an increase both in the oscillation period and in the 
negative peak amplitude after the pulse. The dip in experi- 
mental electric field before the pulse is not reproduced by the 
application of the filters to the simulated data. 

The filters were applied to simulated spectra at carrier trap- 
ping times varying from Tc = 0.1 ps to Tc = 100 ps. As shown 
in Fig.lHl this enables a direct comparison between the exper- 
imental /peak as a function of ion dose, and /peak extracted 
from the model as a function of I/tc. When Tc = 0.1 ps 



the filtered spectral peak is at 2.1 THz, in excellent agree- 
ment with the value (2.2 THz) of the highest dose sample 
(2.5 X lO^^cm-2 at 1 MeV, 1 x lO^'' cm"^ at 2.4MeV). As 
the carrier trapping time is increased (I/tc decreased) the 
simulation tends to /peak ~ 1.6 THz, in agreement with the 
experimental value. Values of the FWHM of filtered simu- 
lated spectra are approximately independent of Tc, owing to 
the bandwidth limitations imposed by the 0.2 mm-thick ZnTe 
crystal. 

A smaller increase in /peak was observed over the ex- 
perimental range of vacancy concentrations: at 10^^ cm~"^ 
/peak = 1.6 THz, and at 10^* cm"^ /peak = 1.9 THz. This 
suggests that carrier trapping may be more significant than 
carrier-vacancy momentum scattering in determining THz 
pulse duration. 

The experimental reduction in the spectral peak of THz 
emission after annealing (Fig.|2J may be attributed to a com- 
bination of a reduction in carrier-vacancy scattering and an 
increase in carrier trapping time from Tc ~ 0.1 ps to Tc 1 ps. 



V. CONCLUSION 

We have studied the surface field THz emission from GaAs 
implanted with multiple doses of high energy arsenic ions both 
experimentally and via a three-dimensional carrier dynamics 
simulation. The effect of trapping defects and neutral impu- 
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FIG. 8: (Colour online) Frequency of peak emitted THz power 
/p (circles, solid line) and FWHM of spectra (triangles, solid 
line) as a function of ion dose at 2.4 MeV (the 1 MeV ion 
dose was a quarter of the 2.4 MeV dose), /p is defined as the 
midpoint of the two frequencies at which the amplitude of 
the Fourier transform of the electric field is 95% of its max- 
imum. The experimental data can be compared with values 
extracted from spectra produced by our carrier dynamics sim- 
ulation (dashed lines), plotted as a function of inverse carrier 
trapping time I/tc. The trend in /p (crosses) from the simu- 
lation reproduces the experiment. The FWHM from the sim- 
ulation (diamonds) do not vary with I/tc, and are consistent 
with those measured. 
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rities on carrier dynamics was included in the model. At a 
short carrier trapping time of 0.1 ps our model predicts an 
increase in the FWHM of THz emission of nearly 80% above 
the FWHM at long (> 10 ps) trapping times. This bandwidth 
increase was experimentally observed as a shift in peak spec- 
tral power towards higher frequencies at greater ion implant 
doses. We anticipate similar bandwidth improvements for 
photoconductive switches made from our GaAs:As^ samples, 
but with greater available THz power. The momentum scat- 
tering rate of carriers from vacancy defects was found to be 
~ 10^'*s~^, comparable to that of the two dominant scattering 
mechanisms in unimplanted GaAs (electron-hole and carrier- 
charged impurity scattering), for the vacancy concentrations 
expected in the highest dose sample. While the THz pulse 
duration decreased at high vacancy concentration, carrier- 
vacancy scattering was not found to account for the full shift 



in peak experimental spectral power. While the simulation 
parameters were chosen to correspond to our GaAsiAs^ sam- 
ples, the model is applicable to LT-GaAs, in which arsenic 
antisites (Asca) are also the dominant defect, and which has 
comparable trapping times. 
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